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THE GRAPH SPECTRUM OF BARYCENTRIC 
REFINEMENTS 

OLIVER KNILL 


Abstract. Given a finite simple graph G, let Gi be its barycentric 
refinement: it is the graph in which the vertices are the complete sub¬ 
graphs of G and in which two such subgraphs are connected, if one is 
contained into the other. If Ao = 0 < Ai < A 2 < • • • < A n are the eigen¬ 
values of the Laplacian of G, define the spectral function F(x) — A[ nx ] 
on the interval [0,1], where [r] is the floor function giving the largest 
integer smaller or equal than r. The graph G 1 is known to be homotopic 
to G with Euler characteristic x(Gi) = x(G) and dim(Gi) > dim(G). 
Let G m be the sequence of barycentric refinements of G = Go- We 
prove that for any finite simple graph G, the spectral functions i 7 b m 
of successive refinements converge for m -A 00 uniformly on compact 
subsets of (0,1) and exponentially fast to a universal limiting eigenvalue 
distribution function Fa which only depends on the clique number re¬ 
spectively the dimension d of the largest complete subgraph of G and 
not on the starting graph G. In the case d — 1, where we deal with 
graphs without triangles, the limiting distribution is the smooth func¬ 
tion F(x) = 4sin 2 (7nr/2). This is related to the Julia set of the quadratic 
map T(z) =4 z — z 2 which has the one dimensional Julia set [0, 4] and F 
satisfies T(F(k/n)) = F(2k/n ) as the Laplacians satisfy such a renor¬ 
malization recursion. The spectral density in the d — 1 case is then 
the arc-sin distribution which is the equilibrium measure on the Julia 
set. In higher dimensions, where the limiting function F still remains 
unidentified, F' appears to have a discrete or singular component. We 
don’t know whether there is an analogue renormalization in d > 2. The 
limiting distribution has relations with the limiting vertex degree distri¬ 
bution and so in 2 dimensions with the graph curvature distribution of 
the refinements G m . 


1. Introduction 

The spectral theory of graphs parallels to a great deal the 

spectral theory of compact Riemannian manifolds USE]. There are areas 
with a good match like inverse spectral topics, heat kernel related topics like 
Hodge theory where the dimension of the space of harmonic solutions Lf = 0 
for example is the number of connected components and more generally, har¬ 
monic fc-forms are related to the k’th cohomology groups in a rather explicit 
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way as applying the heat flow to k forms converges to harmonic forms. This 
can be used for example to construct the Hurewicz homomorphism from 
7 'Tfc(G) to the cohomology groups H k {G ) by applying the heat flow e~ lLk on 
a fc-form / supported initially on a fc-sphere representing an element in the 
homotopy group i r&(G). It leads to a harmonic fc-form representing a homol¬ 
ogy class. Discrete McKean-Singer [26] str(e _tL ) = y(G) illustrates further, 
how the graph theory parallels the continuum m- the super trace of the 
heat kernel for t = 0 is the super trace of 1 which is the definition of the 
Euler characteristic, while in the case t oo, it counts the alternating sum 
of the Betti numbers, which by Euler-Poincare is the same. The cohomology 
of discretizations of manifolds has since the beginning of the development of 
algebraic topology been used to compute the cohomology of the manifold. 
For the convergence of the spectrum of p-form Laplacians of a compact con¬ 
nected Riemannian manifold M of dimension d > 2, Mantuano’s theorem 
[39] tells that for any e discretization G given by a graph, the spectra of 
the Hodge Laplacians L p are related by cX^^iG) < A k,p(M) < CA/^G), 
where the constants c, C only depend on dimension d, the maximal curva¬ 
ture and radius of injectivity of M. Also for eigenfunctions, the number of 
nodal regions of the k eigenfunction /& is bound by a theorem of Fiedler 
m by £;, paralleling the Courant nodal theorem. The ground state energy 
Ai > 0 is estimated in the same way from below by the Cheeger constant. 
The nodal regions on geometric graphs correspond to Chladni figures on 
compact Riemannian manifolds. Also in topology like for Brouwer-Lefshetz 
m or Jordan-Brouwer theory [333, the notions translate nicely to the dis¬ 
crete and barycentric refinements help. There are also places, where things 
are different: the Weyl law relating the growth of the eigenvalues with vol¬ 
ume has no direct discrete analogue because the spectrum of a graph is a 
finite set. Similarly, the Minashisundaram-Plejel zeta function J2\>o^~ S °f 
a manifold needs analytic continuation in the continuum, while in the graph 
case, the zeta function is an entire function. Already in the case of a circle, 
where the Dirac version of that zeta function is the Riemann zeta function, 
the zeros of the discrete analog on circular graphs can be analyzed well, 
unlike in the continuum, where it it is the Riemann hypothesis. The conver¬ 
gence of the roots [36] of the graph zeta functions is a limit when looking at 
barycentric subdivisions of circular graphs, which is also important for Ja¬ 
cobi matrices generalizing graph Laplacians mm which we studied in [22] 
in an ergodic setup, where the limiting operators are almost periodic Jacobi 
matrices over the von Neumann-Kakutani system. The later is the unique 
fixed point of the 2:1 integral extension operation in ergodic theory. The 
hull of the operators is the compact Abelian group of dyadic integers, where 
space translation is addition by 1 and where the renormalization step is the 
shift. We look now at higher dimensional analogue questions to these one 
dimensional Jacobi spectral problems. There are other analogies between 
Schrodinger operators and geometric Laplacians: while isospectral defor¬ 
mations are possible in the former case, geometric Laplacians are harder 
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Figure 1 . Barycentric refinements of the triangle G — K%- 
The number of vertices of the ball G m grows exponentially 
like 0(6 m ) and is exactly known. 


to deform as the isospectral set of geometries is discrete in general. Hav¬ 
ing worked with isospectral deformations of Jacobi matrices in mm we 
looked in m at higher dimensions and noted that writing L = D 2 can en¬ 
able isospectral deformation and break the spectral rigidity. Much later, we 
realized that in the Riemannian as well as in the graph case, one can deform 
the exterior derivative d Eoi m in an isospectral way. The Riemannian 
case is a bit more technical in that story as the deformed exterior deriva¬ 
tives dt are pseudo differential equations in the continuum but they satisfy 
df = 0 so that they deform cohomology. The deformation does not affect 
the Laplacian L, but in the complex, the nonlinear deformation becomes 
asymptotic to the wave evolution. Interestingly, both in the discrete as well 
as in the continuum, space expands, with an inflationary start. This is not 
artificially placed into the model, but is a basic property if one lets the Dirac 
operator move freely in its symmetry group. Graphs can serve as a labo¬ 
ratory to test concepts related to physics. It is reasonable for example to 
see the evolutions of the various differential forms under the wave dynamics 
as manifestations of different forces. It is a caricature for physics, where 
nothing has put into the system except for geometry, which is a Riemannian 
manifold in the continuum and a graph G in the discrete. When answer¬ 
ing the question which graph to take, symmetry again could give a hint: 
take a scale invariant graph. While impossible to achieve in a finite discrete 
structure, taking a graph which is multiple way refined in a barycentric way 
comes close. Not only are parts of the graph homeomorphic to the entire 
graph, but also the nerve graphs of open covers are homeomorphic to the en¬ 
tire graph. This motivates to look at refinements and its spectral properties. 

In pH], we looked at a product of graphs which satisfies the same Kiinneth 
formula for cohomology as in the continuum. This product of G and H is 
defined as follows: take two disjoint unions of G and H. Pick a complete 
subgraph x of G and a complete subgraph y of H. Let xy be the complete 
subgraph generated by rU y. These points xy form the vertices of G x H. 
Two such vertices xy and uv are connected if one is contained in the other. 
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Figure 2. Barycentric refinements of the octahedron. 

A special case is if H = K\ in which case G\ = G x K\ is the barycentric 
refinement of G: its vertices are the complete subgraphs x of G and two 
complete subgraphs x, y are connected if one is contained in an other. The 
graph G\ is homotopic to G. It especially has the same Euler characteristic. 
The dimension of G\ is bounded below by the dimension of G. If G is 
geometric, then G\ is even homeomorphic to G and the dimensions of G and 
G\ are the same. For any d-dimensional geometric G, the graphs G m are 
all (d + l)-colorable, the color being the dimension of the original simplices 
which make up now the vertices. The automorphism group of G also acts 
on G\ so that fixed points of graph automorphisms m can be realized as 
vertices in G\. Barycentric refinements again stress the point of view taken 
by discrete Morse theory nans] that complete subgraphs in a graph can 
be treated as “points”. 

For practical triangulations, barycentric refinements are rarely used. The 
reason is the huge difference between the possible vertex degrees, as each 
refinement doubles the maximal degree. But the appearance of large ver¬ 
tex degrees has also benefits: for any m, the G m graphs are Eulerian with 
Eulerian spheres and have the property that a natural geodesic flow can be 
defined on them in any dimension. This is essential, as spheres and lines are 
important in any geometry and that for defining lines, we need an Eulerian 
property as with an odd degree vertex, the continuation of a “straight line” 
is ambiguous: we have difficulty for example to continue a line through a 
vertex of an icosahedron as we will have to chose what direction to continue. 
For graphs in which all spheres and all spheres in spheres are Eulerian, there 
is a natural continuation. Graph products and especially barycentric refine¬ 
ments have this property. The large and nonuniform degree of some density 
distributed vertices has the advantage as it allows to emulate the rotational 
symmetry in the discrete. For a d-dimensional graph, and considering G m , 
we can in dimension d start with a vertex degree of 0((d + l)! m ) direc¬ 
tions. Unlike in regular lattices or tessellations, the spheres appear more 
rounded and the fact that the unit sphere S(x ) has asymptotically a sim¬ 
ilar amount of directions than the number of vertices in G m gives us the 
property that the exponential map exp x from the unit sphere to the graph, 
covers most of the graph. This is probably the closest we can get to a 
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Figure 3. Barycentric refinements of the tetrahedron G — 
K 4 . The number of vertices grows like 0(24 m ). We know 
V 4 (G m ) = 24 m but we don’t have an explicit formula yet for 
the number of vertices vq(G m) if Go = K 4 . 


Hopf-Rynov property in the discrete without digressing into the quantum 
(as nature does) and look at the wave equation utt — Lu for the Hodge 
Laplacian L — (d + d*) 2 , which because L — D 2 has the d’Alembert solution 
u(t) — cos(Dt)u(0) + sin (Dt)D~ l u r (u), a superposition of wave group solu¬ 
tions e ±lDt = cos (Dt) d= isin(Dt) of the Schroedinger equation ?// = 
with i/j = u + iD~ x v! for the Dirac operator D = d + d* (see |2B]) for which 
Hopf-Rynov is just linear algebra in the graph case as the unitary group is 
finite dimensional and for any two vertices x,y one can solve the problem 
to start a geodesic on e x and reach e y by choosing the correct velocity u'(0) 
and time t. 

Iterated barycentric subdivisions have long been used in topology and ho¬ 
mology theory. The notion is mostly defined in geometric realizations used 
to refine a given triangulation of a simplex in Euclidean space. Note that we 
don’t look at any embedding in an Euclidean space but look at Barycentric 
refinement in a completely combinatorial setup. There are also probabilistic 
connections mm, as iterated barycentric subdivisions of a triangle define a 
random walk on SL(2, R ) and are linked to concrete geometric properties like 
angles triangles. Such work could suggest that unlike in the one dimensional 
case, in higher dimensions, higher dimensional group actions and random 
walks could play a role when studying the spectrum. In graph theory, where 
no embedding into an Euclidean space is required, barycentric subdivisions 
are known in the context of flag complexes or Whitney complexes in partic¬ 
ular. We are not aware of any work on the nature of the spectrum of such 
graphs. 

The inductive dimension of a graph G is defined inductively as the average 
of the dimensions of all unit spheres, incremented by 1 [251 EH- A graph G 
is a geometric d-dimensional graph, if every unit sphere is a (d — l)-sphere 
or a {d — l)-ball and the boundary set of vertices with a ball unit sphere 
form a (d — 1)-dimensional geometric graph without boundary. Geometric 
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graphs play the role of manifolds with or without boundary. The definition 
of graph theoretical spheres and homotopy are both due to Evako. We have 
defined spheres independently in H31CS] but realized in p2| that it is al¬ 
ready in the work of Ivashchenko=Evako [TT1 T6l 0, [TO] . The definition of 
these Evako spheres is based on Ivashchenko homotopy which is homotopy 
notion inspired by Whitehead [48] but defined for graphs. A simplification 
done in [5] makes this homotopy much more practical, like for Lusternik- 
Schnirelman theory for graphs m, to define what a homeomorphism for 
graph is [32] or for graph colorings GUI ESI, where constructing an Eulerian 
3 sphere having a given 2 sphere embedded could explain why 4 colors are 
enough to color a planar graph. 

We proved in [35] the inequality dim(G x H) > dim(G) + dim(iT), which 
holds for all finite simple graphs G, H. This implies in a special case 
dim(Gi) > dimG), still for all finite simple graphs G. In the geometric case, 
the dimension of G\ is the same as the dimension of G. The barycentric 
refinement process G G\ honors geometric graphs, graphs for which the 
unit spheres are spheres. Starting with a d-simplex Go = JQ+ 1 , successive 
refinements G m produce geometric graphs which are d-balls, for which the 
boundary is a (d— l)-sphere. In the case of a triangle G = K% for example, 
the graph G\ is a wheel graph with 6 spikes. It has one interior point and 6 
boundary point. Successive refinements produce larger and larger balls for 
which the boundaries are spheres. 

If G is a geometric d-dimensional graph, then G\ is homeomorphic to G 
in the sense of [32] : there is an open cover Uj of G\ such that the nerve 
graph of Gi is G. To see this, start with unit disks Vj centered at the ver¬ 
tices of G\ which originally were vertices. Then add to Vj the d-simplices 
on the outside which contain (d — 1) simplex in Vj and a vertex simplex 
z. This produces Uj. Two different Uj intersect in a d-simplex, if and only 
the original vertices were connected. Homeomorphic graphs have the same 
topological properties like cohomology, fundamental groups or connectivi¬ 
ties. We furthermore know that for any starting point G, the dimension of 
G m converges to the dimension of the largest complete subgraph in G [35] 
and the reason is the same as for the spectral convergence studied here: the 
largest dimensional simplices dominate. 


2. Counting 

The refinements G m grow exponentially fast already in the 1-dimensional 
case, where the number of edges doubles exactly in each step. 

Lemma 1. a) For G = K 2 , vo(G m )l + 2 m and v\(G m ) — 2 m . 

b) For G = K%, the number of vertices vo (G m ) = 1 — 3(2 m_1 + 2 m + 

2 m ~i 3 m —1) ; the number of edges is Vl (G m ) = 3(-2 m " 1 + 2 m + 2 m " 1 3 m ) 
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and the number of triangles V 2 (G m ) = 6 m . 

c) For G = K d , we have v d (G m ) = ((d+ l)!) m and vq(G m ) > ((d+ l)!) m_1 . 


Proof, a) We have (uo(0), uo(l), uo(2),...) = (2, 3,5,9,17,33,...) as vo(m + 

= 1 . 

'0 2 0 1 

0 2 6 0 


1) = 2vo (m) — 1. Furthermore v\ (m) + vq (rn 
b) From the affine recurrence given by A = 


0 0 6 0 
0 0 0 1 


and S~ l AS = 


Diag(6, 2,1,0) with S 


formula 


1111 
3 10 0 
2 0 0 0 
0 0 10 


we get the exact simplex data 


SE m S~ 


1 

CO 

1 _ 


' 1 - 3 • 2 m_1 + 3 • 2 m + 2 m “ 1 3 m ) ' 

CO 


3(—2 m_1 + 2 m + 2 m-1 3 m ) 

1 


gm 

1 


-1 

T—1 

_1 


which in the first coordinate gives the above sequence 3, 7,25,... which 
grows like 6 m /2. 

(vo(0),v 0 (l),v 0 (2),...) = (3,7,25,121,673,3937,23425,...) • 


c) follows from the fact that in each step, a d-dimensional simplex is cut into 
(d + 1)! simplices. This shows v d (G m ) = (d + l)\v d (G m -i). Furthermore 
^ 0 (Gm) = Vj(Gm- 1 ) > Vd(G m - 1 ). □ 


In |36j, a combinatorial interpretation of the eigenvalues 1, 2, 6 are given. We 
hope it will be possible to get also exact formulas for Vk(G m ) if Go = JQ+1 
if d > 3. To do so, one could use the handshaking relation (see [23] ) 

^2 v k (x) = v k+1 (k + 2 ) , 

x&V 

where Vk(x) are the number of fc-dimensional complete subgraphs 14+1 i n 
the unit sphere S(x). In the case k — 1, this is the Euler handshake giving 
the number v\ of edges in terms of the vertex degrees Vo(x) of the graph. In 
the graph G m , the spheres are all joins of smaller-dimensional spheres and 
in particular come in finitely many types. What one would have to do first 
is to give formulas for Vk(Sd,i,m)) f° r these d-spheres of type l in level m and 
then get ^(G m ) by the Handshaking lemma using the fact that we know 
the number of these spheres. For example, in the case d — 2, only degrees 
4, 6, 8,12,16, 24,.. can occur in the interior and degrees 3, 5, 9,17,... at the 
boundary. It looks like a nice but not hopeless combinatorial challenge to 
analyze this. 
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Some numbers. 

1 ) In the one dimensional case G — where vq(G m ) — 2 m + 1 and 
v\(G m ) = 2 m , we have 

v1 o(G rn+ i)/vo(G rn ) —^2 = 2!. 

2) In the two dimensional case G — we have everything explicit. In 
particular, vq(G m ) = 1 — 32 m_1 + 32 m + 2 m_1 3 m ) which satisfies 

v o(G m +i) / vo(G m ) — 7-6 = 3 !. 

3) In the three dimensional case G — K 4 , we have the following data for 
v = (vo,vi,v 2 ,v 3 ): 

v(G) = (4,6,4, 1 ) 
v{Gi) = (15,50,60,24), 
v(G 2 ) = (149,796,1224,576), 
v(G 3 ) = (2745,17000, 28080,13824), 
v(G 4 ) = (61649,...,.. M 331776). 

We don’t have a closed formula for the sequence vq(G m ) starting with 
4,15,149, 2745,.... We only know that vo(G m -\-i)/vo(G m ) -A 24 = 4!. 

4) In the four dimensional case G = we have computed so far: 
v(G) = (5, 10 , 10 , 5, 1 ), 

u(Gi) = (31,180, 390, 360,120), 

v(G 2 ) = (1081,11340, 33300, 37440,14400). 

v (G 3 ) = (97561,1728000). 

Also here, we only know vo(G m +i)/vo(G m ) —>■ 120 = 5! and have no closed 
formula for the sequence vq(G m ) starting with 5, 31,1081, 97561, ... yet. 

3. Spectral distribution functions 

The Laplacian of a graph is the matrix L = B — A, where A is the adjacency 
matrix and B is the diagonal degree matrix. For a complete graph K^+i for 
example, the matrix is everywhere —1 except for the diagonal, where the 
entries are d. As subtracting d+ 1 produces a matrix B with d dimensional 
kernel and trace — d — 1, the eigenvalues of the Laplacian of the complete 
graph Kd+i is d + 1 with multiplicity d and 0 with multiplicity 1. For the 
house graph G for example, where a triangle K% (roof) is glued on top of a 
square C 4 leading to a graph of dimension 22/12 as it is mixture of one and 
two dimensional components, the Laplacian is 


2 

-1 

0 

-1 

0 

-1 

CO 

-1 

0 

-1 

0 

-1 

CO 

-1 

-1 

-1 

0 

-1 

2 

0 

0 

-1 

-1 

0 

2 


While experimenting with nodal surfaces for eigenfunctions of the Laplacian 
on discrete 3-spheres, where we expect the nodal surface of the ground state 
to be a 2-sphere in general, we also computed eigenvalues of refinements and 
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got interested on how the eigenvalue distribution depends on refinements. 
While we expected some limit to be reached, we would have thought that 
the limit to depend on the topology of the initial graph. This is not the case. 
Starting with a 2-torus or a 2 -sphere produces the same limiting function, as 
the figures illustrate and as we will show below. We can not push numerical 
experiments far, as the number of vertices of the refinements grows so fast. 
But the theory confirms this. The proof is rather elementary using a Lidskii 
lemma. Here is the main result: 

Theorem 2 (Central limit of barycentric refinements). The sequence Fq m (x) 
converges in L 1 ([0,1]) to a limiting distribution function Fd(x ) which only 
depends on the dimension d of the largest complete subgraph of G. 

The convergence is uniform on compact subsets in [0,1) and exponentially 
fast in L 1 ([ 0 , 1 ]) norm or in the L°°([a, b]) norm for every compact interval 
[a, b] C (0,1). In the case d — 1, where we can compute the spectrum 
for circular graphs, we have an explicit limiting function F . We first show 
that the limit FQ m (x) exists for each graph and then prove that the limit 
is universal and only depends on the dimension d of the largest complete 
subgraph of G. The fact that a sequence of monotone functions f m in 
L 1 ([ 0 , 1 ]) which converges in L 1 ([ 0 , 1 ]) converges also uniformly on every 
closed interval [a, b] with 0 <a< 6 <lis known in real analysis (e.g. [58]). 
That the convergence can not be pushed to the boundary point x = 1 follows 
from: 

Lemma 3. For d > 2, the function Fj is not in L°°([0,1]), as the values 
FcmiX) grow exponentially with m. 

Proof. The Courant-Fischer estimate (which is a special case of the Schur 
inquality in linear algebra) shows that \ n -i > max x (deg(:r)), if n is the 
number of vertices. Therefore, for every x: then max(i?,Lr)/(r,r) = L xx = 
d. Since the maximal degree increases indefinitely, the maximal eigenvalue 
X n -i = Fb m ( 1) grows, where n = n(m) is the number of vertices of G m . □ 

We do not know yet about the nature of the limiting functions in the case 
d > 2 : 

Problem: What is the spectral nature of the limiting density of 
states fi = F'? Does it have a discrete or singular part for d >21 

4. The one dimensional case 

Proposition 4. For d — 1, the limiting function is 4 sin 2 ( 7 nr) so that the 
limiting function F is smooth. 

Proof. The eigenvalues of the Laplacian of circular graphs C n are explicitly 
known using discrete Fourier transform conjugating the Laplace operator L 
on ( 2 (^ 71 ) to the diagonal matrix with entries 2 — 2 cos (2nk/n) = 4 sin 2 ( 7 r£;/n) 
so that 


X n = 4sin 2 (7r£;/n) . 
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By definition, we have FQ rn (x) ~^ 4 sin 2 ( 7 nr). 


□ 


Proposition 5. For d — 1, the Laplacian of the circular graphs satisfies 
the renormalization map L(G m +i) = </>(T(L(G m ))) ; where T(x) = Ax — x 2 
is the quadratic map with Julia set [0,4] and (j) restricts the matrix to an 
invariant subspace of half the dimension. 

Proof. This is a direct verification for matrices. In the Fourier picture, 
it becomes a double angle formula of trigonometry. Algebraically, it is an 
identity for Jacobi matrices, the Laplacian of the free particle on the circular 
graph. For G — C 4 for example, we have the Laplacian 


L = 


K = 


2 

0 

-1 

0 

0 

0 

-1 

0 


0 

2 

0 

-1 

0 

0 

0 

-1 


2 

-1 

0 

-1 



-1 

2 

-1 

0 



0 

-1 

2 

-1 



-1 

0 

-1 

2 



C$ and form T(K) - 

= 4 K 

- K 2 

-1 

0 

0 

0 

-1 

0 

0 

-1 

0 

0 

0 

-1 

2 

0 

-1 

0 

0 

0 

0 

2 

0 

-1 

0 

0 

-1 

0 

2 

0 

-1 

0 

0 

-1 

0 

2 

0 

-1 

0 

0 

-1 

0 

2 

0 

0 

0 

0 

-1 

0 

2 


This gives 


We see that that there are two 4-dimensional subspaces of R 8 on which K is 
isomorphic to the old operator L. In the above notation, we have expressed 
this as (f>(K) — L. □ 

Remark: The quadratic map T is associated to the Julia set of the map 
x 2 — 2 which is the ” tip of the tail” parameter 2: = — 2+i-O in the Mandelbrot 
set M. The polynomial T is conjugated to a Tschebychev polynomial, which 
by the way are the polynomials for which the Julia set is an interval [2j. 
Alternatively, in the context of interval maps, the map T is conjugated to 
the Ulam map f{x) — 4x(l — x). The density of the eigenvalues follows 
the arcsin probability distribution f{x) = (x(4 — x))~ 1 / 2 / 71 supported on 
[0,4]. It has with cumulative distribution function ( 2 / 7 r) arcsin(A/T/2). The 
absolutely continuous measure /a = f(x)l[Q^dx is the natural equilibrium 
measure on the Julia set, the measure which maximizes metric entropy and 
equals it to topological entropy log(2). 

5. Convergence 

Intuitively, universal convergence is a consequence of the fact that every sim¬ 
plex in G m spans more simplices producing a self similar pattern. Smaller 
dimensional parts and contributions from boundary get washed away as the 
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refinement progresses, as in dimension d, there are exponentially more points 
than in dimension d — 1. Also in higher dimensions, we expect this to be a 
setup for a renormalization scheme, as the spectrum on part of the graph is 
close to the spectrum of the entire graph and different regions grow in the 
same manner. 


We can not write the limiting F as a fixed point of a renormalization map 
yet. It can not be ruled out yet that the scalar Laplacian of G m +i could 
be related with the form Laplacian of G m . In two dimensions, where only 
1-forms are Fermionic, the spectrum of D 2 is determined by the spectrum 
of the scalar 0-forms and the 2-forms. 


The refinements become more and more self-similar, even so the degrees 
become larger and larger. We will see that the average eigenvalue A (G) = 
Ylk=i A k/ n — 11 Fq 11 1 converges for m -A oo in such a way that A(G m +i) — 
A (G m ) decreases with an exponential rate depending only on the dimension 
of the largest complete subgraph. 

Problem: Is there also in dimension d > la functional equation 
for which the limiting F is a solution? Is there a relation with an 
equilibrium measure of a Julia set as in d = 1? 

When we look at convergence of Fq^ in L 1 ([0,1]), the limiting graph den¬ 
sity 2vi/vo relating the number of edges v\ with the number of vertices vo 
matters. 

Lemma 6. ||Fb m ||i = Fc rn (x) dx = 2vi(G m )/vo(G m ) which is the aver¬ 
age vertex degree of Gm. 


Proof. The trace of the matrix L is the sum of the eigenvalues. By the Euler 
handshaking lemma, it is twice the number v\ (G m ) of edges. The average 
eigenvalue is therefore 

2vi(G m )/vo(G m ) • 


□ 


In the case of a circular graph, we have ||i 7 b m || + 1 = 2. In the case of the 
triangle graph, we have 


II F g , 


3 2 n+1 (3 n + 1) 
3 2 n + 6 n + 2 


-A 6 . 


The fact that the average degree goes to 6 follows of also from Gauss-Bonnet, 
as K(x) — 1 — d(x)/6 is the curvature for graphs without K 4 subgraphs. 
Since the Euler characteristic of the triangle is 1, the sum of the curvature 
has to converge to 1. More generally, we have: 


Corollary 7. In any dimension d, the one has ||F( 3 m ||i -A (d + 1)! expo¬ 
nentially fast. 
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Proof. Indeed, as the number of boundary simplices grows like ( d\) n and the 
number of interior simplices grows like ((d+ l)!) n , the convergence is of the 
order l/(d + l) n . □ 

Since the largest vertex degree grows exponentially in m, we know 11 Fq^ | |oo ^ 
oo. The Cheeger inequality for graphs allows to say something about the 
ground state energy, the smallest nonzero eigenvalue of the connected graph 
G. 

Corollary 8. The ground state energy X\(G m ) converge to 0 exponentially 
fast for rri oo. 

Proof. In dimension 1 , where we know the explicit spectrum and because 
the number of vertices grows exponentially. In general, it follows from the 
fact that the Cheeger number \C(f)\/mm(\A(f)\\B(f)\) goes to zero expo¬ 
nentially fast. □ 

In order to compare eigenvalues, one could use a perturbation result of Weyl 
which tells that for selfadjoint matrices A , £>, the eigenvalues A k(A + B) are 
sandwiched between A k(A) + Ai(£>) and A k(A) + A n (B). Since we need not 
only the individual eigenvalues to converge but also need a l\ -convergence, 
the following variant of the Lidskii’s theorem which was used already in [235] 
is better suited: 

Lemma 9 (Lidskii). For any two selfadjoint complex n x n matrices A,B 
with eigenvalues aq < «2 < • • • < a n and fd\ < /?2 < • • • < fin, one has 
£"=i \ a i ~Pi\< X”j=i \A-B\ij. 

The reduction to the standard Lidksii theorem was told me by Yoram Last 

m- 

Proof. Denote with 7 i £ R the eigenvalues of the selfadjoint matrix C := A— 
B and let U be the unitary matrix diagonalizing C so that Diag( 7 i,..., 7 n ) = 
UCU*. We calculate 

Em = = ^{-l^UikCkiUu 

i i i,k,l 

< Ei^i-iEv^^e/I^Ei^i- 

k,l i k,l 

The claim follows now from Lidskii’s inequality JY |ctj — (3j\ < \ jj\ (see 

m) □ 

Proposition 10. For any d > 1 and G = JQ+ 1 , there exists a limiting 
eigenvalue distribution Fd(x) = lim m ^ 00 Fq^ (x). 

Proof. We proceed by induction. For d — 1, we have explicit eigenvalues. 
For d — 2, we make a refinement and divide it up into 6 pieces. Lidskii shows 
that the eigenvalues of G m+ i consist of 6 copies of the eigenvalues of G m plus 
a correction term which comes from lower dimensional interfaces and is much 
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Figure 4. For the proof, we cut the d-dimensional graph 
G m into (d+ 1)! smaller pieces of the form G m -The mod¬ 
ification on the boundary affects exponentially less vertices 
so that by Lidskii, the l\ difference between the spectra goes 
to zero exponentially fast with m. This applies then to the 
P([ 0 , 1 ]) functions F Gm . 


smaller. In words, ||F Gm+12 || = ||U ; m , 2 II + 6 • 3||F Gm l ||(2 n /6 n ). As ||F Grri l || 
converges, the sequence ||Tb m2 || is a Cauchy sequence. Now, lets go to the 

case d = 3. We have |T Gm+lj3 || = H-Fg^sII + 24 '3|TG m , 2 ||(6 n /6 n2 ) and 
again have a Cauchy sequence showing that the limit exists in PQO, 1]). In 
each dimension we get a limiting function as the lower dimensional interfaces 
between the similar chambers grow exponentially slower than the chambers 
themselves. □ 

Corollary 11. There exists a limiting density of states /id = F' d (x) which 
is a measure on [0,1]. 

Proof. As Fcm are monotone, it defines a measure /i m . Since the limiting 
eigenvalue distribution is monotone, it defines a measure fa. □ 

With pointwise convergence, we would have weak convergence of /i m to (i . 
But we don’t know that yet. 

Lets now look at the proof of the theorem: we have to show that the limiting 
function F does not depend on the initial graph. Having seen convergence 
for a simplex K&+ 1 , it follows for a finite union of k simplices glued along 
lower dimensional simplices. Each simplex evolves in the same way under 
the barycentric evolution and the average is the same function F. When 
cutting the graph apart or disregarding lower dimensional parts, the modifi¬ 
cations lead to change on a set of vertices which becomes exponentially less 
relevant as m grows as it is lower dimensional. This is justified by Lidskii’s 
estimate as we can estimate the sum of the eigenvalue differences. 
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6. Vertex degree distribution 

Lets look at the vertex degree distribution d\ < d 2 < ds .... In the same 
way as for the spectrum, we can define a degree distribution function 

HGm ( X ) d[ xn \ . 

The degree and eigenvalue distributions are linked by the Schur’s inequality 
so that one can deduce the following from the eigenvalues. It is also possible 
to do it directly: 

Corollary 12 (Limiting vertex distribution). There exists a limiting vertex 
distribution function H(x) — lim m ^ oc (x) . It is independent of the 
initial graph G and depends only on the dimension d of the largest complete 
subgraph of G. The distribution satisfies H{x) > F(x) forx G [0,1], H( 0) = 
F(0),H(1) = F(1). 

Proof. Splitting the graph G m into (d+ 1 )! smaller graphs G m - 1 and noting 
that the degrees of the lower dimensional walls can be neglected in the limit 
as the number of vertices in those parts grows slower, we have a Cauchy 
sequence in L 1 . 

An inequality of Schur For any selfadjoint matrix A, the Schur inequality 

t t 

E a * - E ^ 

2=1 2=1 

holds, where fii are the eigenvalues of A and cq are the diagonal elements 
a*i > ai+ 1 of A. In the case t — n we have equality since tr(L) = JT X The 
Schur inequality gives Yll di > Yll =i ^i if fh e degrees d{ of the graph and 
eigenvalues \{ of the graph are ordered in an ascending way. □ 

Remark. A bit stronger than Schur is a result of Gone |1] ) which assures 
that 

k k 

di < Mu 1 ^ k < n — 1 
2=1 2=1 

if Mi > /i 2 > • • • > v n = 0 are the eigenvalues of the Laplacian of G and 
di > d 2 > • • • > d n > 0 are the vertex degrees of the graph G. This im¬ 
plies that (x) > Fq m (x) for each m and suggest that also in the limit 
H(x) > F(x) for every x G (0,1). 

In the case d — 2, the degree d{ is related to the curvature K(x) = 1 — d(x)/6 
which adds up to the Euler characteristic of G. We see that the curvature 
distribution is related to the eigenvalue distribution. 


In general we can look at the curvature 


k(x) = ,-YM + YM-YM 


2 


3 


4 
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Figure 5. To the left is a graph of the spectral density func¬ 
tion of the Dirac operator d + d* of a barycentric refinement 
of the house graph. To the right an illustration of the Schur 
inequality: the integrated vertex degree function is an upper 
bound for the integrated eigenvalue function F(t) dt and 
agrees at the end points 0,1. 


of a vertex which by Gauss-Bonnet-Chern [24] adds up to the Euler charac¬ 
teristic x(G) of G. 

The limiting distribution works also for the Dirac operator D — d + d* and 
so for the Hodge Laplacian D 2 = dd* + d*d. The argument is the same. 

Corollary 13. For any k, the Laplacian L^(G m ) on k-forms has a spectral 
limiting function which only depends on k and the dimension d of the largest 
simplex in G = Go. 

In each dimension d, the graph G m belonging to a simplex Go = dQ can 
be cut (using some modifications on smaller dimensional sub graphs) into 
(d+1)! smaller isomorphic pieces, which are of the form G m -\. The (d— 1)- 
dimensional cuts do not matter in the limit as the number of simplices in 
them grows exponentially less fast. 

Since the matrices D 2 (G m ) and L(G m + 1 ) have the same size, we compared 
the distribution of the Hodge Laplacian on the level m with the distribution 
of the scalar Laplacian on the level m + 1 and they appear comparable. Be¬ 
cause of super symmetry, the 1-form Laplacian L\ has the same spectrum 
than the union of the 0 and 2-form Laplacians. The 2-form Laplacian is 
related to the Laplacian of the dual graph. 

There might be more relations if things are extended to Schrodinger op¬ 
erators. In m , we looked at renormalization maps L -A D satisfying 
L — D 2 + c, where D is a Laplacian on a barycentric refined circular graph. 
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Its not clear whether there are algebraic relations for d > 2. For a trian- 

2 -1 -1 “ 

-1 2 -1 
-1 -1 2 

refinement G\ and the form Laplacian D 2 (Gq) have the same size 


gle G = Gq with Laplacian 


the Laplacian L{G\) of the 


L(G0 = 


D {Go) = 


CO 


-1 

0 

-1 

-1 

0 
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-] 
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7. Figures 



Figure 6 . The spectrum of the barycentric refinements 
G m of the triangle, for m = 3,4, 5. 
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Figure 7. The spectral functions F in the case d = 1, d = 2 
and d — 3. We see the spectrum of barycentric refinements 
Gj for G = C 4 , the spectrum of the refinement G 4 of the 
octahedron G and finally the spectrum the refinement G 3 of 
the tetrahedron K 4 . 



Figure 8 . Barycentric refinements of the house graph. The 
one dimensional component can be neglected in the limit. 



Figure 9. The spectrum of barycentric refinements G m of 
the house graphs for m = 3,4, 5, which is the case d — 2. 
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Figure 10. The spectrum of barycentric refinements G 4 
for the triangle the octahedron and the torus. The eigen¬ 
value distribution look almost identical. These are all 2 di¬ 
mensional cases. 




Figure 11 . A refinement of the triangle and its dual graph: 
the graph whose vertices are triangles and where two trian¬ 
gles are connected if they intersect in an edge. The dual 
graph does not have triangles. 



Figure 12. The Hodge Laplacian L — (d + d *) 2 for refine¬ 
ments of the triangle K% has three blocks Lo,Ti,T 2 , where 
Lk acts on fc-forms. The Hodge Laplacian of the tetrahedron 
K 4 has blocks Lq, Ti, £ 2 ? ^ 3 - Super-symmetry assures that 
the union of the spectra of the Bosonic parts Lo,T 2 is the 
union of the spectra of the Fermionic parts Li, L 3 . 
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